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O 1 ABSTRACT 

Aims. We study the molecular environment of the luminous blue variable star G24. 73+0.69 to 
investigate the origin of the two infrared shells around this massive star and determine its effects 
CNj ' on the surrounding interstellar medium. 

Methods. We analyze the distribution of the molecular gas using the 13 CO J= 1—0 emission ex- 



> 

traded from the Galactic Ring Survey. We use near and mid-infrared data from 2MASS and 
GLIMPSE to identify the young stellar objects in the field. 

Results. We discover the molecular counterpart to the outer infrared shell around G24.73+0.69. 
The CO shell was probably blown by the stellar wind of the star mainly during its main sequence 
phase. We also find molecular gas that corresponds to the inner infrared shell, although its origin 
remains uncertain. We identify seven young stellar objects within the molecular material, whose 
birth might have been triggered by the stellar wind of the luminous blue variable star. We sug- 
gest that both G24.73+0.69 and the progenitor of the nearby supernova remnant G24.7+0.6 were 
formed from the same natal cloud and represent the most evolved members of a so far undetected 
cluster of massive stars. 

Key words, stars: individual: G24.73+0.69 - stars: winds, outflows - ISM: clouds - stars: for- 
mation 



1. Introduction 

Luminous blue variable (LBV) stars are very massive objects that evolve from O-type main se- 
quence (MS) stars burning hydrogen in their core to become Wolf-Rayet (WR) helium core burn- 
ing stars (Maeder 1989; Humphreys et al. 1989). This transitional stage is characterized by a high 
mass-loss rate (typically between 10~ 5 and 10~ 4 M /yr) sometimes accompanied by so-called giant 
eruptions, and significant photometric variability on timescales from months to years (Humphreys 
& Davidson 1994). The LBV phase is short-lived and as a consequence only 12 cases and 23 candi- 
dates have been reported in our Galaxy so far (Clark et al. 2005). As a result of the mass loss, most 
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of the LBV stars (either confirmed or candidates) are surrounded by a nebula that expands with 
characteristic velocities of between 30 km s~'and 200 km s _1 (see Clark et al. 2005). These nebulae 
have dusty and gaseous components (of a few solar masses) and a large variety of morphologies, 
from circular to bipolar (Nota et al. 1995; Weis 2001). 

Since nebulae around LBV stars are strong emitters in the infrared (IR) regime, most of them 
have been observed in this wavelength range. At present, there are only a few cases in which the 
chemical evolution and the presence of circumstellar molecular gas related to these nebulae have 
been determined by molecular studies, the most representative sources being AG Car (Nota et al. 
2002) and G79. 29+0.46 (Rizzo et al. 2008; Jimenez-Esteban et al. 2010). Molecular material was 
also discovered around the yellow hypergiant (YHG) stars IRC+10420 and AFGL2343* (Castro- 
Carrizo et al. 2007). Theoretical models of stellar evolution and observations of post-MS massive 
stars (Meynet & Maeder 2000; Humphreys & Davidson 1994) suggest that low luminous LBVs 
may pass through a cooler red supergiant (RSG)/yellow hypergiant phase prior to the LBV pe- 
riod. The RSG/YHG phase is also characterized by mass-loss activity (for a review, see de Jager 
1998) that can lead to the presence of molecular gas around these massive stars. By analyzing the 
molecular material, we can determine the mass-loss history of the massive star and its effects on 
the surrounding interstellar medium (ISM). 



2. Presentation of G24.73+0.69 

The ring nebula around the LBV star G24.73+0.69 (hereafter G24) was discovered by Clark et al. 
(2003) in the mid-IR. Fig. 1 displays the emission in the 5/?/fzer-IRAC 8 fim band toward G24, 
which clearly shows the nebula around the central massive LBV star. It has a roughly circular 
morphology with a slight elongation in a direction parallel to the Galactic plane. Its shape can be 
approximated by an ellipse of 40"x34". Taking an expansion velocity of 200 km s , Clark et al. 
(2003) estimated a mass-loss rate of ~ 9.5 x lO~ 5 M /yr for the LBV star and an age of ~ 4, 800 
yr for the nebula. They quoted an upper limit of 5.2 kpc for the distance based on the consistency 
between the modeled and the observed spectral energy distribution of the star. In addition, there 
is a fragmentary outer shell (OS) of a bipolar morphology with lobes (labeled 1 and 2 in Fig. 1) 
projecting into opposite directions. This structure forms an incomplete shell of 4.6'x2.3' oriented 
~ 60° with respect to the Galactic plane, whose geometric center does not match the position of 
the LBV star. Fig. 1 also shows the position of the IR dark cloud (IRDC) 024.789+0.633, which is 
adjacent to the lobe 2. All IRDCs are considered as the precursors of protostars, being a signature 
of active star formation (Rathborne et al. 2006, 2007). 

To explain the origin of the IR shells, Clark et al. (2003) argued that the inner one (i.e. the ring 
nebula) originates in material ejected from the central star during the LBV phase. However, based 
on the size and age of the bipolar OS, they suggest that it formed from the interaction between the 
stellar wind and the ISM and/or the natal molecular cloud (MC) in a period prior to the current 
phase. 

In this work, we study the molecular gas toward G24.73+0.69 and its environment with the aim 
of unveiling the origin of these IR shells. 
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Fig. 1. Emission in the Spitzer-JRAC 8 yum band toward the LBV star G24.73+0.69. We indicate the positions 
of the LBV star and its nebula, the lobes of the bipolar outer shell, and the IRDC 024.789+0.633. 

3. Results and discussion 

3.1. The distribution of the molecular gas 

We study the molecular gas around G24 using data from the Galactic Ring Survey (GRS, Jackson 
et al. 2006). The GRS was performed by the Boston University and the Five College Radio 
Astronomy Observatory (FCRAO). The survey maps the Galactic ring in the 13 CO J=l-0 line 
with an angular and spectral resolutions of 46" and 0.2 km s , respectively. The observations 
were performed in both position-switching and on-the-fiy mapping modes, achieving an angular 
sampling of 22". 

From the inspection of the I 3 CO data cube in the whole velocity range, we found morphological 
signatures of a possible association between G24 and the surrounding molecular material in the 
velocity range between +39 and +44 km s 1 (all velocities here being referred to the local standard 
of rest). In Fig. 2, we show the integrated velocity channel maps of the B CO J=l— emission every 
0.7 km s with the mid-IR emission from the Spitzer-IRAC 8 /mi band. From this figure, we can 
discern several condensations of molecular gas forming a molecular shell around G24 (see mainly 
panels at +40.9 and +41 .5 km s _1 ). This shell delineates the exterior border of the bipolar OS. This 
distribution indicates that there is a connection between the molecular gas and the bipolar OS. If 
this were the case, both of them would be at the same distance. Using the rotational model of Fich 
et al. (1989) and adopting a velocity of +42 km s _1 as a systemic velocity for the molecular shell, 
we obtained kinematic distances of either 3.5 kpc or 12 kpc. Lacking any additional discriminator 
between the two kinematic distances and taking into account that a distance of 5.2 kpc is quoted as 
an upper limit for G24 (Clark et al. 2003), we adopted a distance of 3.5 kpc as the most plausible 
value for the LBV star and the molecular shell. This new distance led to a substantial reduction 
in the calculated intrinsic luminosity of G24. Clark et al. (2003) quoted Zog(L*/L ) = 5.6 for a 
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Fig. 2. Two-color image of the field around the LBV star. In red: Spitzer-IRAC 8 jim band. In green: the 
integrated 13 CO J=l— emission every 0.7 km s -1 . The CO contours levels are 0.5, 0.9, 1.8, and 2.7 K km s -1 . 

distance of 5.2 kpc. Taking a distance of 3.5 kpc, the luminosity drops to /og(L+/L ) ~ 5.25, 
which would make it one of the faintest LBVs yet identified (see for comparison the HR diagram 
of Clark et al. 2009), thus a possible post RSG/YHG object. 

We were able to derive a number of parameters characterizing the CO. By assuming local 
thermodynamic equilibrium (LTE), we estimated the mass of the molecular shell using the 
equations 

N( 13 CO) = 2.42xl0 14 1 ^ e ! 5 ° 29 ^ J r n dv, 

where T ex is the excitation temperature of the 13 CO transition and T13 is the optical depth of the 
line. Assuming that the 13 CO J=l-0 line is optically thin, we can use the approximation 

/T]T.d\> ~ f Tndv, 

J(T ex ) - J(T h ) J 

where 

5.29 

^ ' ~ e 5.29/T _ l ' 

Th = 2.7 K is the background temperature, and Tb is the brightness temperature of the line. We 
assumed that T ex = 10 K and used the relation N(H 2 )/N( 13 CO)~ 5 x 10 5 (Simon et al. 2001) to 
estimate the column density of molecular hydrogen N(H2). The mass of the molecular clumps was 
calculated from 

M^m H ^[D 2 nN(H 2 )], 

where D. is the solid angle subtended by the 13 CO J=l-0 beam size, m u is the hydrogen mass, fi is 
the mean molecular weight assumed to be 2.8 by taking into account a relative helium abundance 
of 25 %, and D ~ 3.5 kpc is the distance. Our summation was performed over the area of each 
molecular clump. We summed the mass of the clumps around G24 and found that the molecular 
shell has a total mass of ~ 1, 000 M . The original ambient density no, an important parameter to 
study the dynamics of the shell, could be estimated by assuming that the calculated total mass was 
originally distributed in the volume limited by the bipolar OS, which has a mean radius of ~ 3.5'. 
At a distance of 3.5 kpc, this corresponds to a mean radius ro ~ 3.5 pc. Using this value, we found 
that «o ~ 230 cm -3 . 
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3.2. The origin of the molecular shell 

There is compelling evidence of shells of interstellar material produced by the effect of the strong 
stellar winds of massive stars on the ISM. These so-called stellar wind bubbles have been observed 
in neutral (Cappa et al. 2005; Cichowolski et al. 2008; Giacani et al. 201 1) and molecular (Cappa 
et al. 2009, 2010) gas. We explored the hypothesis that the action of the stellar wind of the LBV star 
and/or its progenitor O-type star has blown a cavity of interstellar material piling up the molecular 
gas in the border and forming the molecular shell discovered in the present work. 

To test this scenario, we compared the mechanical energy released by the star into the ISM 
during the lifetime of the molecular shell, with the kinetic energy Ek of the swept-up material. We 
first estimated the dynamical age (?</„,) of the shell using the equation of Whitworth et al. (1994) 
for a stellar wind bubble: td yn ~ Q.Qlr^L^^n 1 ^ Myr, where ro is the radius of the shell in pc, 
L37 is the mechanical luminosity L w of the stellar wind in units of 10 37 erg s _1 , and 113 is the 
original ambient density in units of 10 3 cm -3 . We calculated L w = l/2Mv^, by taking the typical 
parameters of an O-type star in Prinja et al. (1990) and Mokiem et al. (2007), namely a mass-loss 
rate M — 2 x 10~ 6 M Q /yr and a stellar wind velocity v w - 2,000 km s _1 . Thus, we found that 
L w ~ 2.5 x 10 36 erg s _1 . Taking the mean radius of the shell ~ 3.5 pc and the original ambient 
density ~ 230 crrT 3 , we evaluated the td v „ to be ~ 0.2 Myr. From this result, we can see that the 
CO shell is more than an order of magnitude older than the age of the LBV phase, which typically 
has a duration of < 10 4 yr. This indicates that the detected CO shell may be the material swept-up 
by the stellar wind of the central massive stars during its MS phase. 

The kinetic energy can be estimated as = l/2MVg, where M is the mass of the shell and 
v e its expansion velocity. We calculated the expansion velocity to be v e ~ 6.2L^ 5 n 3 " 1, ' 5 f^ 5 km 
s -1 (Whitworth et al. 1994), where t dyn = 0.2 Myr. We obtained v e ~ 13 km s _1 . Using M = 
1, 000 M Q , we inferred that £^ ~ 1.7 X 10 48 erg. For the mechanical energy, we assume that E l0 = 
Lojtdyn ~ 1- 6xl0 49 erg. We calculated the efficiency of the energy conversion of the wind e = ^ and 
obtained a value of ~ 0.11. This agrees with the expected value for an energy-conserving stellar 
wind bubble (Mellema & Lundqvist 2002), hence we conclude that the central star has blown a 
stellar wind strong enough to create the molecular shell. We also evaluated the contribution of the 
star during the LBV phase to the energy outflow. Taking M = 9.5 x 10~ 5 M Q /yr and v w = 200 km 
s _1 from Clark et al. (2003) and ~ 10 4 yr as a typical duration for the LBV phase, we obtained 
E^f v ~ 3.8 x 10 47 erg, which represents less than 3 % of the energy released by the star during the 
MS phase. 

These calculations link the formation of the molecular shell to the stellar wind of the central star 
during the MS phase, in good agreement with the origin of the bipolar OS suggested by Clark et al. 
(2003). In this scenario, the elongated morphology of the bipolar OS (and also of the molecular 
shell) may be a consequence of the asymmetric stellar wind from the central star. Evidence of 
non-spherical winds has been found in other massive stars. For example, Meaburn et al. (2004) 
discovered a giant lobe projecting from the ring nebula around the LBV star P Cygni. According 
to the authors, the lobe formed from the mass-ejection activity prior to the current eruptions that 
produced the nebula. They suggested that the peculiar morphology of this giant lobe is due to the 
interaction between the ejected material and an asymmetric cavity in the ISM blown by the toroidal 
stellar wind of the central star during the MS phase. 



A. Petriella et al.: The molecular gas around G24.73+0.69 

Finally, we note the presence of molecular gas superimposed on the ring nebula (see Fig. 2, 
panels at +42.1 and +42.8 km s _1 ). This material may i) originate in a mass-ejection event from 
the central LBV star, as observed in the ring nebula around G79. 29+0.46 (Jimenez-Esteban et al. 
2010); or ii) be part of the molecular shell. Using the equations of Sect. 3.1, we estimated a mass 
of ~ 20 M and a density of ~ 120 cm 4 . The obtained mass is higher than the mass of the nebula 
around G24 estimated by Clark et al. (2003) (~ 0.45 M ) and of the nebulae around other LBVs 
(see Smith & Owocki 2006 and Clark et al. 2009). Thus, we propose that the molecular material 
located right upon the ring nebula belongs mainly to the molecular shell. However, we cannot rule 
out the presence of some molecular gas ejected by the LBV star. We plan to observe a field around 
G24 with higher density tracers and higher resolution to confirm the presence of this material. 

3.3. Star formation activity around G24. 73+0.69 

The environment of the LBV star G24 might be one that promotes triggered star formation. As 
we have shown in the previous section, the massive star has been interacting with the neighboring 
molecular gas by means of its strong stellar wind. To search for primary tracers of star formation, 
we identified the young stellar object (YSO) candidates using the GLIMPSE Point Source Catalog 
(GPSC) in the Spitzer-IRAC bands. We constructed a [5.8]-[8.0] versus [3.6]-[4.5] color-color (CC) 
diagram for the sources around G24. In the CC diagram, we used the photometric criteria of Allen 
et al. (2004) to classify the sources according to their evolutionary stage: class I are protostars with 
prominent circumstellar envelopes, class II are disk-dominated objects, and class III are the most 
evolved YSOs whose emission comes mainly from the central star. 




24.850 24.800 24.750 24.700 24.650 
Galactic longitude 



Fig. 3. Two-color image toward G24.73+0.69. In green with contours: emission of the 13 CO integrated be- 
tween +39 and +44 km s~'(contour levels are 2.7, 4.4, and 6.2 K km s~'). In red: 8 /jm band. The magenta 
crosses are the class I YSO candidates selected according to the criterion discussed in the text. 

Fig. 3 displays the distribution of the class I YSO candidates around G24. We can discern seven 
of these YSO candidates superimposed on the molecular gas swept-up by the massive star. Four of 
them (sources 1, 2, 3, and 4) spatially coincide with the maximum of molecular emission. We note 
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that this region coincides with the IRDC 024.789+0.633 (indicative of active star-forming activity) 
and lies next to the lobe 2 of the bipolar OS. 

We fit the spectral energy distribution (SED) of the class I YSO candidates using the fluxes in 
the 2MASS J, H, and Ks bands and in the four Spitzer-IRAC bands. We used the tool developed by 
Robitaille et al. (2007) and available online 1 . We assumed the distance to be between 3 and 4 kpc. In 
Table 1, we report the class I YSO candidates magnitudes in the 2MASS and Spitzer-IRAC bands. 
Robitaille et al. (2006) defined three evolutionary stages based on the values of the central source 
mass M+, the disk mass Mdisk, the envelope mass M env , and the envelope accretion rate M env . Stage 
I YSOs are those that have M env /Mi, > 10 yr~ l , i.e., protostars with large accretion envelopes; 
stage II are those with M^/M* > 10~ 6 yr~ x and M^/M* < 10~ 6 yr~ l , i.e., young objects with 
prominent disks; and stage III are those with M^/M* < 10~ 6 yr~ l and M env IM* < 10~ 6 yr~ l , i.e., 
evolved sources where the flux is dominated by the central source. The goodness of the fitting for 
each model was measured by ax 2 value. We defined the "selected models" as those satisfying the 
equation^- 2 - xt- < 2Af, where xt- corresponds to the best-fitting model and N is the number of 
input data fluxes (fluxes specified as upper limit do not contribute to AO- The fitting tool also fit the 
IR fluxes to a stellar photosphere to check whether the source could be modeled by a main-sequence 
star with interstellar reddening. The goodness of the fitting is defined by ax\ ar value. 

Table 1. 2MASS and Spitzer-IRAC magnitudes of class I YSO candidates. 



Source 


Source (GPSC) 


2MASS Qual. 


/ 


H 


Ks 


3.6 //m 


4.5 jum 


5.8 /im 


8.0 pm 


1 


G024.7956+00.6558 


N/D 








13.047 


12.173 


11.715 


11.321 


2 


G024.805 1+00.6456 


UUA 


18.131 


17.059 


13.707 


10.548 


9.589 


8.791 


8.059 


3 


G024.7901+00.6410 


N/D 








12.874 


11.096 


9.983 


9.639 


4 


G024.7907+00.6379 


N/D 








13.439 


12.518 


11.638 


11.059 


5 


G024.7545+00.6771 


AAA 


14.818 


13.227 


12.158 


10.967 


10.497 


10.194 


9.859 


6 


G024.7532+00.6400 


UUA 


14.643 


12.784 


12.312 


11.792 


11.463 


11.722 


11.386 


7 


G024.7735+00.6068 


AAA 


13.898 


13.172 


12.803 


11.907 


11.485 


10.083 


8.577 



Notes. 2MASS Qual.: A is the best photometric quality, with a snr > 10, and U means that the magnitude 
value is an upper limit. N/D indicates no detection of the source in the 2MASS survey. 



In Table 2, we report the results of the SED fitting for the seven class I candidates within the 
molecular shell. For all of them, all the selected models correspond to stage I and II (except source 
2 for which we get one stage III model). Moreover, for all of them apart from source 6 we get 
Xmin < x\ar- Thus, we can confirm that they are young sources at the earlier stages of evolution. 

As we note in the previous section, the molecular shell might have formed as a consequence of 
the interaction between the stellar wind of the central massive star and the ISM. The presence of 
YSO candidates within the molecular shell and an IRDC are signs of active star-forming activity. 
We suggest that the birth of these new stars might have been triggered by the expanding wind 
bubble acting on the molecular gas. 



1 http://caravan. astro, wise .edu/protostars/ 
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Table 2. Results of the SED fitting for class I YSO candidates. In the last column, we report the evolutionary 
stage of the source, based on the criterion of Robitaille et al. (2006) for the selected models. When the models 
correspond to different stages, the less likely is presented in brackets. 



Source 


X 2 ■ IN 

71 mm' 


xlJN 


Stage 


1 


0.4 


15 


I or (II) 


2 


4 


77 


II 


3 


0.1 


93 


I 


4 


0.02 


37 


I or (II) 


5 


0.9 


4 


II or (I) 


6 


2.3 


2 


II 


7 


10 


91 


II 



3.4. G24.73+0.69 and its surrounding: the big picture 

Massive stars form in giant molecular clouds which, under certain circumstances, fragment into 
dense clumps that collapse to form new stars. These stars spend a short period of time (~ 10 6 yrs) 
in the MS before evolving to become LBV and WR stars, which eventually end their lifetimes 
exploding as supernovae. These explosions usually occur when the progenitor star still lies close 
to its natal cloud and while other companion stars may still be in the MS or LBV/WR phase. As 
a consequence, several supernova remnants have been observed interacting with evolved massive 
stars (see Velazquez et al. 2003; Sushch et al. 201 1) and molecular gas (see Jiang et al. 2010 for an 
exhaustive list of SNR/MC associations). 

The LBV star G24.73+0.69 lies ~ 7' from the SNR G24.7+0.6. In the radio band, this remnant 
displays a couple of incomplete shells and a polarized filled central core with a flat spectrum (Reich 
et al. 1984), which indicates that it consists of a plerion powered by an undetected pulsar. Thus, 
the SNR formed from the collapse of a massive star. Petriella et al. (2008) reported the existence 
of molecular gas interacting with G24.7+0.6 in the velocity range between +38 and +50 km s _1 , 
which interestingly places the SNR at the same kinematic distance as G24. In addition, active star 
formation probably triggered by the SNR and/or its progenitor was also discovered around the 
remnant (Petriella et al. 2010). In Fig. 4, we show the emission of the 13 CO J=l-0 integrated in 
the velocity range between +38 and +50 km s . We indicate the position of the different features 
in the field: the LBV star (red cross), the SNR (in blue contours), and the MC interacting with the 
remnant (cloud 2 in the nomenclature of Petriella et al. 2008). 

From the previous figure, we see that this large region, which harbors a rich variety of objects 
(a LBV star and its nebula, a SNR, abundant molecular material hosting dark clouds, possible star- 
forming sites), is an ideal laboratory to investigate the genetic connection between different stellar 
populations. On the basis of the morphology and distribution of the molecular gas, it is likely 
that the progenitor star of the SNR G24.7+0.6 and the LBV star may have formed from the same 
giant MC that now we observe fragmented into smaller features. They may be the most evolved 
members of a so far undetected cluster/association of massive stars that also formed from the same 
molecular material. In addition, they may have triggered the formation of a second generation of 
stars that now we observe as young objects deeply embedded in the molecular gas. This scenario 
should be investigated further by searching for the missing intermediate-age stellar population. 
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Fig. 4. Emission of the 13 CO integrated between +38 and +50 km s . The blue contours are the radiocontin- 
uum emission at 20 cm of the SNR G24.7+0.6 taken from the MAGPIS. The red cross indicates the position 
of the LBV star G24.73+0.69. 

namely massive stars still in the MS phase, which formed together with the LBV and the supernova 
progenitor. 

4. Summary 

We have analyzed the 13 CO emission in the surroundings of the LBV star G24.73+0.69. We have 
discovered a fragmented molecular shell in the velocity range between +39 km s _1 and +44 km s _1 , 
which delineates the infrared bipolar outer shell of G24. On the basis of spatial, morphological, and 
kinematic coincidence of features, we suggest that there is a connection between this IR shell and 
the molecular gas. We argue that the molecular shell formed from the interstellar material swept-up 
by the stellar wind of the central star mainly during its MS phase. The elongated morphology of 
the bipolar and molecular shells may be a consequence of the asymmetric stellar wind from the 
central massive star. We have also detected molecular emission probably associated with the inner 
infrared nebula, but observations of higher angular resolution are needed to spatially resolve this 
emission and establish its origin. 

We have studied the star formation activity around G24 and discovered seven YSOs with spec- 
tral characteristics of protostars, projected on the molecular cloud that interestingly coincides with 
the IRDC 024.789+0.633. Thus, this region harbors the typical components commonly found in 
the vicinity of star forming regions, namely an IRDC, several solar masses of molecular gas, and 
embedded YSO candidates. We propose that the birth of these young objects might have been 
triggered by the expanding stellar wind bubble. 

From the study of the distribution of the molecular gas in a large region toward G24, we suggest 
a link between the origin of the LBV star and the progenitor star of the nearby SNR G24.7+0.6. 
This SNR is also interacting with the neighboring molecular gas and shows star-forming activity 
in its surroundings. We propose a scenario where both massive stars belong to a first generation 
of stars that formed from the same natal MC, pointing to the possibility of finding other massive 
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stars with a common origin. The strong stellar winds of these massive stars might have triggered 
the formation of a second generation of stars, which now appear as protostars deeply embedded in 
the molecular gas. 
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